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ABSTRACT 



Results are presented from a detailed analysis of optical and X-ray observa- 
tions of moderate-redsliift galaxy clusters from the Canadian Network for Ob- 
servational Cosmology (CNOC) subsample of the EMSS. The combination of 
extensive optical and deep X-ray observations of these clusters make them ideal 
candidates for multiwavelength mass comparison studies. X-ray surface bright- 
ness profiles of 14 clusters with 0.17 < z < 0.55 are constructed from Chandra 
observations and fit to single and double jS models. Spatially resolved tem- 
perature analysis is performed, indicating that five of the clusters in this sample 
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exhibit temperature gradients within their inner 60 — 200 kpc. Integrated spectra 
extracted within R2500 provide temperature, abundance, and luminosity informa- 
tion. Under assumptions of hydrostatic equihbrium and spherical symmetry, we 
derive gas and total masses within R2500 and R2oo- We find an average gas mass 
fraction of /gas(R2oo) = 0.136±0.004 Yi'^''^ , resulting in VLm = 0.28±0.01 (formal 
error). We also derive dynamical masses for these clusters to R2oo- We find no 
systematic bias between X-ray and dynamical methods across the sample, with 
an average Mdyn/Mx-ray = 0.97 ± 0.05. We also compare X-ray masses to weak 
lensing mass estimates of a subset of our sample, resulting in a weighted average 
of Miens/Mx-ray of 0.99 ± 0.07. Wc investigate X-ray scaling relationships and 
find powerlaw slopes which are slightly steeper than the predictions of self-similar 
models, with an E{z)~^ Lx-Tx slope of 2.4 ± 0.2 and an E{z) M25oo-^x slope 
of 1.7 ± 0.1. Relationships between red-sequence optical richness (Bgc^red) and 
global cluster X-ray properties {Tx, Lx and M2500) are also examined and fitted. 

Subject headings: cosmology: observations — X-rays: galaxies: clusters — galax- 
ies: clusters: general 



Introduction 



Clusters of galaxies are valuable cosmological probes. As the largest gravitationally 
bound objects in the universe, they play a key role i n the tracing and modeling of large scale 



structure formation and evolution (e.g., IVoitI |2004 iBahcall. Fan, fc CenI 119971 ). Because of 



our rapidly expanding sample of known clusters, finding efficient means of estimating cluster 
properties is a highly desirable goal. 

The contribution of cluster studies to the field of observational cosmology hinges on our 
ability to accurately estimate cluster masses. In particular, through the determination of 
both gas mass and total mass, cluster analysis can lead to estimations of the cosmological 
mass density, Qm, while accurate measurement of the evolution of the cluster mass function 
provides important constraints on both the nor malization of the matter power spectrum, 



(78 an d w, the dark energy equation of state (e.g., iLevine. Schulz. fc White! l2002l : lEke et al. 



19981 ). In addition, gaining an understanding of the evolution of X-ray scaling relationships 
(such as Lx-Tx and M-Tx) with redshift provides an important contribution to our abili ty 



to accurately model the evolution of large-scale structure in the universe (e.g., IVoitll2004l ). 



Intermediate redshift (0.2 < z < 0.6) clusters are well situated for the study of cluster 
properties. They are compact enough to be observed without tedious mosaicing, and they are 
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present in statistically significant numbers. Intermediate redshift clusters are also luminous 
enough to permit detailed investigation, with the potential for placing strong contraints on 
cluster evolution. 

Three frequently applied approaches to estimating cluster masses are gravitational lens 
modeling, optical spectroscopic measurements of the cluster galaxy velocity dispersion, and 
characterizing X-ray emission as a means of tracing the underlying potential well of the 
cluster. Each of these methods uses different observations and assumptions, which can 
be tested through their direct comparison. An alternate optical approach to the efficient 
estimation of cluster masses involves the use of optical richness. This method relies on the 
assumption that galaxy light is a reliable tracer of total cluster mass, and requires calibration 
via other methods. 

Lensing mass estimates test both the assumption of hydrostatic equilibrium and our 
knowledge of the mass distribution in clusters. Because they probe all of the projected 
mass along the line of sight, which may inclu de additional mass concentrations, they ar e 
susceptible to overestimation of cluster mass (ICenI 119971 : iMetzler. White, fc LokenI I2OOII ). 
Dynamical mass estimates work under the assumption that velocity dispersion is directly 
related to the underlying gravitational potential of the cluster. Pitfalls of this technique also 
include the danger of overestimating the clu ster mass in cases where substructu re or mergers 
drive up velocity dispersion measurements (IBird. Mushotzky. fc Metzlerlll995l ). 



The hot, diffuse intracluster medium (ICM), which is observed in the X-ray, should be a 
direct tracer of a cluster's underlying potential well. Under assumptions of isothermality and 
hydrostatic equilibrium, the surface brightness of a cluster of galaxies can provide information 
on gas density as well as total gravitating mass. Factors that influence the accuracy of X-ray 
mass estimates are temperature gradient s, substructure, and rnergers , which can compromise 
the previously stated assumptions (e.g.. iBalland fc BlanchardI 119971 ). 



One of the main objectives of this work is to identify, through direct comparison, any 
systematic biases in these methods with the ultimate goal of determining a robust calibra- 
tion between optical richness and cluster mass. Optical richness measurements are easily 
available due to the fact that their estimation requires very little observing time. Therefore 
an accurate calibration of the relationship between mass and optical richness would allow 
us to determine the masses of large samples of clusters in a highly efficient manner, provid- 
ing strong constraints on the evolution of the cluster mass function and, consequently, key 
cosmological parameters. 

In this paper w e present a detailed analysis of h igh resolution Chandra X-ray observa- 
tions of 14 CNOC (lYee. EUingson. &: Carlberglll996l ) clusters at z ~ 0.3. In Sections 
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we probe the temperature, metallicity, morphology, and surface brightness of the hot ICM 
present in each cluster. From this initial analysis, we derive mass estimates (Section [5]) 
which are compared to dynamical and weak lensing mass estimates (Section [7]), and then we 
examine the X-ray scaling laws of our sample (Section [H]). Finally, in Section [HI we use our 
results to calibrate relationships between red-sequence optical richness (Bgc,red) and global 
cluster X-ray properties (Tx, Lx and M25oo)- 

Unless otherwise noted, this paper assumes a cosmology of Hq = 70 km s~^ Mpc~^, 
Qrn = 0.3, and Qa = 0.7 



2. Cluster Sample & Chandra Observations 



X-ray observations of our sample were retrieved from the Chandra Data Archive (CDA) 
after conducting a search for current ly available Chandra observations of the Canadian Net- 
work for Observational Cosmology (lYee. EUingson. fc CarlbereJEooel. C NOC) intermediate 
redsh ift (0.17 < z < 0.55) sub sample of 1 5 EMSS clusters (ICioia et al. Ill990l : iHenry et al. 



1992 ) an d one Abell cluster (Abell 1958h. The selection criterion for this sample can be 
found in Yee. EUingson. fc Carlberg (jl996 ). This sample has been extensively observed by 
the CNOC cluster survey (CNOC-1), and galaxy redshifts of ~ 1200 cluster members as well 



as detailed photometric catalogues are available for these clusters (e.g.. lYee et al. Ill996l ) 



Chandra Advanced CCD Imaging Spectrometer (ACIS) observations of 14 of the 16 
CNOC clusters were obtained. Six of these clusters were observed with the ACIS-S CCD 
array, and eight were observed with ACIS-I, with an overall range in individual exposures 
of 11.7 — 91.9 kiloseconds. Three of the clusters in our sample were observed on multiple 
occasions, in which case we have chosen the longest of those observations to include in our 
analysis. Each of the observations analyzed in this study possesses a start date which falls 
on or after 24 April, 2000, indicating a focal plane temperature of —120" C for the entire 
sample. 

Aspect solutions were examined for irregularities and none were found. Background 
contamination due to charged particle flares were reduced by removing time intervals during 
which the background rate exceeded the average background rate by more than 20%. The 
event files were filtered on standard grades and bad pixels were removed. A two-dimensional 
elliptical Lorentzian was fit to the counts image of each dataset to locate the center of 
the X-ray emission peak. All centroid position errors are within a resolution element (~ 
0.5"). In the case of the heavily substructured cluster MS0451+02, a fitted central peak was 
unobtainable, so a position at the center of the extended emission was chosen for spectral 
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analysis at an RA, Dec of 04:54:09. 941, +02:55: 14.52 (J2000), and the surface brightness 
profile was centered on the slightly offset peak of emission at an RA, Dec of 04:54:07.249, 
+02:54:27.31 (J2000). 

Table [1] provides a list of each of the clusters in our sample, including redshifts, obs-id, 
detector array, and corrected exposure information for each observation. 

After initial cleaning of the data, 0.6-7.0 keV images, instrument maps, and exposure 
maps were created using the CIAO 3.2 tools DMCOPY, MKINSTMAP and MKEXPMAR 
Data with energies below 0.6 keV and above 7.0 keV were excluded due to uncertainties in 
the ACIS calibration and background contamination, respectively. 

Flux images were created by dividing the resulting images by their respective expo- 
sure maps. Point source detection was performed by running the tools WTRANSFORM 
and WRECON on the flux images. Adaptively smoothed flux images were c reated with 



CSM QQTH. Figured] contains CFHT images obtained from the CNOC database (lYee. EUingson. fc Carlbei 



19961 ). and, where available, HST images obtained from the MAST archive, each of which is 



overlayed with X-ray contours created from smoothed Chandra flux images. 



3. Surface Brightness 

Using the 0.6-7.0 keV images and exposure maps, a radial surface brightness profile was 
computed in 1" annular bins for each cluster. These profiles were then fit with both single 
and double (3 models. Single P models take the form 



I(r) 



/b + /o 1 + -7 



where Is is a constant representing the surface brightness contribution of the background, 
Iq is the normalization and rc is the core radius. 



The double /9-model has the form 
I(r) =I^ + lJl + 



-3/3i 



1 + 



-3/32 + ^ 



(2) 



where each component has fit parameters {In,rn, Pn)- 



Though a fair number of the clusters in our sample posess detectably elliptical or irreg- 
ular emission, circula r surface brightness pro f iles w ere c hosen based on a numbe r of factors. 
First, as discussed in lNeumann fc Bohringen (119971 ) and iBohringer et al. I (119981 ). differences 
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in masses derived using elliptical vs. circular (3 profiles are small (~ 5%). Secondly, cir- 
cular j3 profiles provide a more straightforward comparison between dynamical and X-ray 
derived masses, since dynamical mass calculations assume sphericity. Finally, temperature 
uncertainties should outweigh any errors introduced by assuming radial symmetry. 

The parameters of the best fitting single (3 model of each cluster are shown in Table 
|2l While nine of the clusters' surface brightness profiles do not exhibit excess unmodeled 
emission in their cores when fit with a single (3 model, five of the clusters do exhibit excess 
emission in their cores. The surface brightness profiles of these clusters are better fit with 
the addition of a second (3 component. All of the clusters requiring a double j3 model are 



believed to c ontain cooling flow s : Abe 



MS1455-F22 flLewis et al 



1999 



Allen 



.1 2390, MS0440+02, MS0839+29, MS1358+62, and 
20001) and exhibit strong (LHr>- )- j N n] > 10"^^ erg s"^) 



extended Ha emission ( iDonahue. Stocke. &: Gioial Il992l : iLewis et al. I Il999l ) . Best fitting 



double (3 model parameters are given in Table [31 Surface brightness profiles of the clusters 
in our sample are shown in Figure [2l Discrepancies in background values are due to the 
use of both ACIS-I and ACIS-S arrays (ACIS-S typically has background values which are 
a factor of three higher than ACIS-I). 



4. Spectral Analysis of X-ray Observations 

4.1. Integrated Spectra and R2500 

For the purpose of fitting cluster scaling relationships within a well defined region of mass 
overdensity, global cluster properties were determined for our sample within R25oo- Using the 
fitted X-ray centroid positions obtained in Section [2] (with the exception of MS0451+02, as 
noted in Section [2]) , a spectrum was extracted from each cleaned event file in a circul ar region 



with a 300 h^Q^ kpc radius. These spectra were then analyzed with XSPEC 11.3.2 (lArnaud 



19961 ). using weighted response matrices (RMFs) and effective area files (ARFs) generated 
with the CIAO tool ACISSPEC and the latest version of the Chandra calibration database 
(CALDB 3.0.1). Background spectra were extracted from the aimpoint chip as far away 
from the cluster as possible. Each spectrum was grouped to contain at least 20 counts per 
energy bin. 

Spectra were fitted with single temperature spectral models, inclusive of foreground 
absorp tion. Each spectrurn was initially fit with the absorbing column frozen at its measured 



value (jPickey fc LockmanI Il990l ). and redshifts were fixed throughout the analysis. Metal 
abundances were allowed to vary. Data with energies below 0.6 keV and above 7.0 keV were 
again excluded. 
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In s ome cases evidence of excess photo electric absorption has been seen in coohng flow 
clusters (lAUen fc FabianI 119971 : lAUenl |2000| ) . To investigate this possibility a second fit was 
performed, allowing the absorbing column to vary. No evidence for excess absorption was 
found in 13 of the 14 clusters. The MS0440+02 spectrum, however, was equally well fit by 
a model inclusive of excess foreground absorption. 

The results of these fits, combined with the best fitting f3 model parameters from Sec- 
tion [31 were then used to estimate the value of R2500 for each cluster. This is accomplished 
by combining the equation for total gravitating mass 



Mtnf (< r) 



T{r)r f 6 In p ^ 6 In T 



GprUp \5 \ni (5 In r 



(3) 



(where pnip is the mean mass per particle) with the definition of mass overdensity 



(4) 



where the critical density pc = SH'^/StcG, z is the cluster redshift, and A is the factor by 
which the density at ta exceeds the critical density. These equations are then combined with 
the density profile implied from the (3 model (assuming hydrostatic equilibrium, spherical 
symmetry, and isothermality) 



Pgas{r) = po 



r 

1 + ^ 



-3/3/2 



(5) 



resulting in the equation 



3pT 



Gpmp{A/3)7Tpc{l + zfr^A 



(6) 



(]Ettorill200d : lEttori et al. \\2QQ4 \. 



After the initial estimation of R25005 additional spectra were extracted from within that 
radius, and spectral fitting was performed again. This iterative process was continued until 
fitted temperatures and calculated values of R2500 were consistent with extraction radii. The 
results of this process, including values of R2500! temperatures and abundances are shown in 
Table HI along with 90% confidence ranges. 
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An additional process was carried out for the five cooling core clusters in our sample. 
Using spatially resolved temperature profiles (Section I4.2p . the radius within which cooling 
became prominent in each cluster (Rcutofr) was estimated. A single annular spectrum was 
then extracted for each cluster from that radius to R2500 and fitted with single temperature 
spectral models. The observation of MS0440+02 did not possess enough signal to constrain 
a fit to the cooling flow excised spectrum, therefore it was discarded. The results of these 
fits are shown in Tabled along with inner extraction radii and calculated values of R25oo- 

Overall, our inte grated temperatures are comparable to the previous ROSAT results 
of iLewis et al. I (119991 ). though they are better constrained through Chandra ob servations. 



More recent a nalyses of subsets o: 




ASCA fMatsumoto et al. 



these clusters have been conducted with ROSAT (Mohr et a. 



2000 



Arabadiis. Bautz. fc Garmire 



Ettori et al. 



2002 



Henrvll2004f). and Chandra (lAUen. Schmidt, fc Fabian 



Donahue et al. 2003; Ettori fc Lombardi 12008 



2OO4J ). Our fitted temperatures are consistent within errors of the vast ma- 



jority of the values reported in these papers. Noticeable discrepancies arise occasionally in 
the case of cooling flow clusters (where some of these authors have not excised the cooling 
flow contribution to integrated temperature) and MS1008-1224, for which we consistently 
derive a lower temperature than other studies. This is not surprising considering the irregu- 
lar morphology of this cluster, and may be due to differences in the regions for which spectra 
were extracted. Metallicities for this sample are consistent with an overall value of ~ 0.3 
solar, as is seen in lower redshift clusters. This suggests that ICM enrichment occurs early in 



a clus ter's history, evidence which is corroborated by the work of iMushotzky &: Loewenstein 
fll997h . 



4.2. Spatially Resolved Spectra 

With the goal of elucidating the radial dependence of temperature for the clusters in our 
sample, spectra were extracted within circular annuli spanning the central 150 — 600 h^Q^ kpc 
of each cluster, depending on the signal-to-noise ratio of each observation. The extraction 
regions were sized to include at least 2500 counts per spectrum in the 0.6-7.0 keV band for 
each data set. This number was chosen as a minimum for achieving reasonable temperature 
constraints on the emission in each annulus. 

Extractions were performed after the removal of previously detected point sources (Sec- 
tion [2]), using the CIAO tool ACISSPEC, which incorporates the changing response over 
large areas of the detector. Three of the 14 clusters in this sample did not possess enough 
signal for more than one extraction region, and were removed from the sample during the 
subsequent spatial analysis. Spectra from the eleven remaining clusters were then grouped 
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to contain 20 counts per bin, and were analyzed using XSPEC. Single temperature spectral 
models with fixed galactic absorption and varying abundances were fitted to the data. The 
results of these fits were used to create a temperature profile for each cluster. 

Figure [3] illustrates the eleven resulting temperature profiles. It is clear from these 
figures that significant cooling cores are present in five of the clusters in this sample: Abell 
2390, MS0440+02, MS0839+29, MS1358+62, and MS1455+22. These are the same five 
clusters which required double (3 model fits to their surface brightness profiles (Section [3]). 
The remaining six temperature profiles are consistent with isothermality. 



4.3. The Universal Temperature Profile 



The spatially resolved temperature profiles of the five cooling core clusters were scaled 
by R2500 and T2500 i n an attempt to check thei r cons istency with the universal temperature 
profile proposed by lAUen. Schmidt, fc FabianI (120011 ). The combined data from these five 
clusters were then fit using the functional form T(R )/T25nn = To + Ti[(x/xc)''/(l + (x/xc)^)], 
with X = R/R2500 flAUen. Schmidt, fc Fabia^boOlh . This fit resulted in Tq = 0.397 ± 0.007, 
Ti = 0.56 ± 0.02, Xc = 0.0865 ± 0.0002, and r/ = 1.3 ± 0.2, with a reduced of 1.3. 



These values are all consistent with the best fit found by lAUen. Schmidt, fc Fabian 



( I2OOII ). The resulting function, however, does not asymptotically approach 1 at large radii, 
ostensibly due to the fact that these temperature profiles do not extend all the way to 
R2500- Figure m shows the individual temperature profiles as well as the binned and averaged 
temperature profile of the five clusters, with the best fitting functions overlayed. 



5. Mass Determinations 

Using the results of the spectral fits and (3 model fits, along with Equations [3l [5], and [6l 
gas masses and total masses were calculated out to both R2500 and i?2oo for the clusters in our 
sample. Central densities were determined via an expression relating the observable cluster 
X-ray luminosity to gas density, using emission measures obtained during spectral fitting in 
XSPEC. 

For the nine non-cooling core clusters, single temperature spectra fits and single j3 
model parameters were used. For the five clusters which exhibit significant cooling, spectral 
parameters from cooling core excised spectral fits were used, and (3 and were taken from 
the results of double P model fits. 
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To determine the effects of inner temperature gradients on total cluster mass, least- 
squares fitting was performed on the temperature profiles obtained in Section 14.21 for the 
five cooling core clusters. The resulting parameterizations were included in Equation [31 
and masses were calculated out to the edge of the cooling region (Tabled]). Masses were 
also calculated without the inclusion of this parameterization, and the results compared. 
According to the outcome of this exercise, the inclusion of temperature gradients in mass 
calculations of the five cooling core clusters in our sample would result in, at most, a 0.3% 
correction to the total mass within the cooling region. This correction is negligible compared 
to the other uncertainties in mass calculations and is therefore not included in the final 
results. 

An additional uncertainty is present in mass estimations of MS045H-0250 due to its 
significant irregularity. Choosing a centroid at the center of the extended emission rather 
than one at the most central peak of emission produces a total mass which is greater by 
25%. The beta parameter which results from using this centroid, however, is unusually high 
(1.9). X-ray mass determinations are presented in Tables [6] and [7] along with 68% confidence 
intervals. 



6. Gas Mass Fractions and ^l.^. 

Under the assumption that clusters provide a fair representation of the universe, gas 
mass fractions, fgas, defined as the ratio of cluster gas mass to total gravitating mass, were 
calculated for our sample within R200 and are hsted in Table [3 fgas can be used to calculate 
the cosmological mass density, Qm, via the relation 



fgas (l + O.W-5)^ 



(7) 



( Allen et al. 2002h. where O.lQ/t"'^ represents the baryonic contribut ion from optically lu- 



minous matter (IWhite et al. Ill993 



0.0223;°:°°°^ is adopted (WMAP, ISpergel et al. 
average fgas 
mass density of Q 



Fukugita. Hogan. fc Peebles! Il998l ). A value of Q,b h'^ of 
20071 ). and here we take h = 0.71. Using an 



0.136 ± O.OOihjQ^^ within R200 for our sample, we calculate a cosmological 
0.28 ± 0.01 (68% confidence, with error bars re presenting the forma l 
error). This value is in good agreement with WMAP three year results (ISpergel et al. 1120071 ). 
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7. Mass Comparisons 

7.1. Dynamical Mass Comparisons 

Comparisons b etween X-ray and d ynamical masses for the CNOC sample were previ- 
ously undertaken by lLewis et al. I (119991 ) using ROSAT observations. While they found good 
agreement between the two methods of mass estimation, accurate surface brightness model- 
ing and detailed investigations of cluster temperature gradients were unavailable due to the 
comparitively poor spatial resolution of ROSAT (particularly at moderate redshift). Here 
we utilize the 0.5" spatital resolution of Chandra to improve the accuracy of these initial 
comparisons. 



Detaile d dynamical studies of CNOC cl usters were perfo r med b v lCarlberg et al. I (1l996l ). lBorgani et al 



ploying an improved interloper-removal algorithm, and I Van der Marel et al. I (120001 ) inves- 



( 1l999l ). and IVan der Marel et al. I (j2000f ). ICarlberg et al. I (119961) provides velocity disper- 
sions obtained from CFHT spectroscopy, iBorgani et al. 1 Jl999h adjusts these v alues by em- 



tigates the isotropicity and galaxy distribution of a composite CN OC cluster created via 
dim ensionless scaling. He r e we will primarily draw from the work of iBorgani et al. I (119991 ) 



and IVan der Marel et al. I (120001 ) for our mass estimates. 



Dynamical masses can be calculated from velocity dispersions via the Jeans equation 



M(r) 



(5 In (t3 5 In uii) 

+ . , ^ ' + 213 



5 In r 



(5 In r 



where ar represents the radial velocity dispersion, u{r) is the galaxy number density profile 



and j3 represents the anisotropy of the system. According to IVan der Marel et al. 
the CNOC clusters can be treated as isotropic (i.e. /3 = 0), and i/(r) takes the form 



mm . 



i/[r] 



Voir /a) [1 + r/a] 



7-3 



(9) 



where the length scale of the mass distribution is set by the parameter a, and 7 represents 
the logarithmic power-law slope near the center. The best fitting values of a and 7 for 
isotropicity are 0.224 and 0.75, respectively. This rel ationship was obtained by creating 



a composite CNOC cluster via dimensionless scaling ( Van der Marel et al. Il2000l ). as was 



a plot of a vs. r. This figure (jVan der Marel et al. 



2000l . Figure 2) indicates th a t a is 



not a strong function of r adius out to -R200; and in keeping with both iLewis et al. I (119991 ) 



and ICarlberg et al. I (119961 ) we assume that o"(r) = a. 



Using Equation [HI dynamical masses were calculated for our sample out to i?2oo (as de- 
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termin ed by X-ray parameters in Section H?T]) . using the velocity dispersions of iBorgani et al. 



( 119991 ) ■ These masses were then compared to X-ray derived masses from the previous section 
(Section^]). A weighted average gives an overall dynamical to X-ray mass ratio of 0.97±0.05, 
where the error bar indicates the uncertainty in the mean. Table E] lists both dynamical and 
X-ray derived masses along with their ratio and 68% confidence intervals. Figure [5] is a plot 
indicating the dynamical to X-ray mass ratio of each cluster in the sample. 

The high dynamical to X-ray mass ratios of MS1006.0+1202 and MS1008. 1-1224 may 
be due to overestimated velocity dispersions of these objects, as they both have markedly 
irregular emission (Figure [1]). However we also see evidence of good agreement between 
X-ray and dynamical masses in irregular objects (MS0451. 5+0250), as well as disagreement 
in some regular objects (MS0839. 8+2938), therefore a clear pattern does not make itself evi- 
dent. Likewise, cooling core objects show no obvious systematic departures from consistency. 
Overall, dynamical and X-ray mass estimations for this sample show remarkable agreement. 



Our resulting ratio of dynamical to X-ray masses is consistent with that quoted by lLewis et al 



(119991 ) in their ROSAT study, however the scatter about the mean of our distribution is ~ 5% 



smaller. This decrease in scatter is indi c ative of the improved spatial and spectral resolution 



of Chandra. In addition, iLewis et al. I (119991 ) systematically overestimate the core radii of 



cooling flow clusters, and therefore their masses, another result which is likely due to the 
poorer spatial resolution of ROSAT. 



7.2. Weak Lensing Mass Comparisons 



Weak lensing mass estimates were obtained for sev en of the cluste rs in our sample using 
deep optical observations at the CFHT 3.6m telescope (jHoekstrall2006l . in preparation). The 
model-independent projected mass estimates that we employ in this paper were calculated 
for the inner 500 kpc of each cluster, using a cosmology of Hq = 100, Qm = 0.3 and 
= 0.7. 

To compare X-ray derived masses to weak lensing masses, we calculated a cylindrical 
X-ray mass within 500 kpc. This was done using previously determined f3 models, the 
adopted cosmology (above), and a cyhndrical mass projection out to 10 Mpc from the cluster 
core. A comparison of weak lensing masses to X-ray derived masses can be found in Table [HI 
a plot of mass ratios is given in Figure O Contributions to lensing signal from structures 
along the line of sight may result in masses which are biased somew hat high. Simil arly, large 
scale structure along the line of sight results in increased scatter (lHoekstrall200ll ). Despite 
these possible challenges, the weighted average of our mass ratios gives a weak lensing to X- 
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ray mass ratio of 0.99 ± 0.08, with a reduced of 0.93. Though the distribution in Figure [6] 
appears assymmetric, it is not statistically significant. 



8. X-ray Scaling Laws 

8.1. The Lx - Tx Relationship 

Unabsorbed bolometric X-ray luminosities within R2500 were obtained for our sample in 
the following manner. Unabsorbed X-ray luminosities for the 2-10 keV band were calculated 
in XSPEC during spectral fitting (Section 14.11) . In the case of the five clusters which exhibit 
significant cooling, corrected, integrated luminosities were obtained by fixing X-ray temper- 
atures at the values determined for cooling core corrected spectra (Table To convert 
2-10 keV luminosities to bolometric X-ray luminosities, correction factors were obtained via 
NASA's Portable, Interactive Multi-Mission Simulator (PIMMS) for each individual cluster, 
using a thermal bremsstrahlung model and the spectrally determined temperature of the 
cluster. The resulting bolometric luminosities are given in Table [TOl 

The Lx — Tx d ata were fit using the biva riate correlated errors and intrinsic scatter 



(BCES) estimator of lAkritas fc Bershadyl (119961 ). This estimator allows for r neasurement er- 



(2001 


) and 


Yee and EUinesonl f 


2003 


)• 


to 


Allen. Schmidt. & Fabian 


(2001 


)^ 


ivas 



^2500 



10^4 erg s-i E(z) 



Ci (kT25oo) 



C2 



(10) 



,1/2 



The best fitting values w ere C 



1999 



1 3+°°^ 



where E(z 

and Cy, = 2.4 ± 0.2. W hile this slope is lower th an some estimates (e.g. lArnaud fc Evrard 



Ettori et al. II2OO4I ). it is consistent with both lAUen. Schmidt. &: Fabianl (|2001l ) and lEttori. De Grandi 



2OO2I ). Because it is steeper than the expected self-similar slope of 2 (jVoitll2004l ). it indicates 
modest negative evolution in the Lx — Tx relationship at z ~ 0.3. 



Ettori et al. I (120041 ) examine scaling laws in clusters at moderate to high redshift (0.4 < 
z < 1.3), so here we adopt their definition of scatter for comparison 

{ay = Sj=iAr(^0(7 Lj ohs ~ log Lj^p^cdY/N ). Defined this way, our scatter in luminosity 
ay = 0.20, is smaller than that of Ettori et al. {ay = 0.35); however they have included 
significantly higher redshift clusters in their fit. A plot of the data with our relationship 
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overlaid is presented in Figure [71 



8.2. The Mass-Tx Relationship 



A BCES fit was also pe rformed between cl uster temperatures and mass estimates. This 
relationship takes the form (lEttori et al. II2004J ): 



{E{z)Mtot) = Ci (kT) 



C2 



(11) 



We again use cluster properties which were determined within R25oo- The best fitting values 



were Ci = l.SlnQXlO^^ and = 1.7 ±0.1. This slope is consistent with lAUen. Schmidt. &: Fabian 



fl200lh . lEttori. De Grandi. fc Molendil fl2002h . and lEttori et al. I fl2004f ). and is again steeper 
than t he expected slope of 1.5 f Voitll2004 ). Our scatter in mass is also lower than Ettori et al. 
(I2OO4J ) (cjy = 0.07 compared to cxy = 0.15). A plot of the data with our relationship overlaid 
is presented in Figure [HI 



9. X-ray Properties and Optical Richness 

Optical richness is in essence a measurement of the galaxy excess in the direction of 
a cluster above a certain magnitude limit and within a specific aperture. The particular 
optical richness parameter that is u sed in this work, Bnr., is d efined as the amplitude of the 
galaxy-cluster correlation function (ILongair fc Seldneiill979l ): 



Dgc I 



(12) 



Galaxies are observed as projections on the sky, therefore what is measured is the 
angular two-point correlation function of galaxies, uj{6), where 6 is the angle on the sky. 
uj{6) is approximated by a powerlaw of the form: 



(13) 



(jPavis &: Peebleslll983l : lYee and Lopez-Cruzlll999l ). where Agg is the galaxy-galaxy angular 
correlation amplitude, ue is taken as the the distribution of galaxies around the the center 
of the cluster, and the amplitude is then relabeled as Age- This amplitude can be measured 
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from an image by counting the background-corrected excess of galaxies within a certain 9 to 
a particular magnitude limit. 

Bgf. is then calculated th rough a deprojection analysis which assumes spherical symme- 



try (ILongair fc Seldnerlll979l ): 



B 



gc 



DT-3A„ 



■I^$[M(mo,z)] 



(14) 



where N^g represents the background galaxy counts to apparent magnitude mo, D is the angu- 
lar diameter distance to the cluster redshift z, is an integration constant, and <l>[M(mo, z)] 
is the integrated luminosity function of galaxies to the absolute magnitude M which corre- 
sponds to mo at z. 

One of the challenges in the calculation of Bg^ involves the lack of a complete knowledge 
of the galaxy luminosity function at high redshifts. Thi s uncertainty can be minimized by 
employing the parameter B geared (IGladders fc Yed l2005l ). which is calculated using galax- 
ies in the red-sequence. This parameter is expected to provide a more robust indication 
of c luster mass due to the w ell understood passive evolution of red-sequence cluster galax- 



ies (Ivan Dokkum et al.lll998l ) as opposed to the more unpredictable nature of star forming 



populations. Since one of our main goals in this work is to provide a comparison sample for 
future studies of clusters at higher redshift, we use the parameter Bgc^red for the measurement 
of richness. 

At the relatively low redshifts of the CNOC clusters, the difference between Bg^ and 
Bgcred is Small. We estimate Bgc^red values by applying sma ll corrections to the origi- 



nal B „r. values of the CNOC clusters (lYee and EUingsonl |2003| ) , based on their blue frac- 



tions (lEUingson et al. Il200ll ). Most of the corrections are of the order of ~ 10%. Values 



of 5, 



gc,red 



(in units of h^^Q 



Mpc 



1.77\ 



as well as the cluster X-ray luminosities obtained in 



Section 14.11 are given in Table [101 To keep the same scale as previous work using the B 
parameter, they are computed using Ho = 50 km s~^ Mpc~^. 



Bgc has previou sly been shown to correla te strongly with the X-ray properties of clus- 



ters in this sample (lYee and EUingsonl l2003l ) . Here we re-examine these correlations us- 
in g improved X-ray data from Chandra. Extending the simple relationships expressed 



m 



Yee and EUingsonl (120031 ) to Bgc^red we have 



B 



gc,red 



OCT 



'7/2 



(15) 



and using our best fit Lx —Tx and Mass— Tx relationships (Section [8]), in combination with 
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Equation [15] 

Bgcred OC ( -^TTK ) ' (^6) 



T \ 7/4.1 



E{z) 



and 



Bg^^red OC (E(^) M2500)^^'-', (17) 

In this section we derive the best fitting relationships between cluster X-ray properties and 
Bgcred for our sample. A generic form for these relationships of 



log X = Ci + C2 log Bgcred, 



is adopted, where X represents the particular property being fit. For T y, E(z) M25nn, and 
L x/E (z), units of keV, M^r^, an d lO'^'^erg s~^ were used, resp ectively. Asin lYee and EUingson 



(l2003l ). the BCES estimator of lAkritas fc Bershadyl (119961 ) was employed. 



The results of these fits are all consistent within errors with those of lYee and EUingson 

1.77). The data, with 



(120031 ) ■ and are also consistent with the e xpected value of 7 (7 



the results of these fits as well as those of lYee and EUingsoru (120031 ) overlaid, are shown in 
Figures [glfTTl Best fitting parameters of all fits are given in Table [TTl 



10. Summary and Discussion 

We have presented a comprehensive analysis of Chandra observations of 14 medium 
redshift (0.1709 < z < 0.5466) clusters of galaxies from the CNOC subsample of the EMSS 
(Table [1]). Imaging analysis has provided information on the relative quiescence of each 
cluster (Figured]). The spatial resolution of Chandra has allowed us to determine the surface 
brightness profiles of these clusters down to 1" scales (Figure[2]). This has enabled us to obtain 
precise models of cluster emission (Table [2]), except in the case of MS0451+02 which displays 
a considerable amount of substructure. Five clusters in our sample (Abell 2390, MS0440+02, 
MS0839+29, MS1358+62, and MS1455+22) exhibit excess emission in their cores, indicating 
the presence of cool gas. The surface brightnesses of these five clusters are better fit with 
the inclusion of a second (3 component to model the excess core emission (Table EJ). 

Chandra's 0.5" spatial resolution has enabled us to study the radial distribution of tem- 
peratures on small scales (~ 7 — 100 kpc) in the inner 150-600 kpc of eleven of these objects 



-17- 



(Figure[3]). While nine clusters in our sample have temperature profiles which are consistent 
with isothermality, the five cooling core clusters (listed above) show clear temperature gra- 
dients in their innermost ~ 60 — 200 kpc. The temper ature profiles of these five a 



consistent with the "universal" temperature profile of I Allen. Schmidt. &: FabianI (120011 ) . 



usters are 



The energy resolution of the instrument has provided well constrained spectral analyses 
of the central R2500 (~ 350 — 760) kpc of each object. Temperatures obtained initially 
from 300 kpc radius extraction regions began an iterative spectral extraction and fitting 
process culminating in the determination of robust integrated temperatures, luminosities, 
and abundances within R2500 (Tables H] and [5]). Cooling core corrected spectra were used 
to determine global temperatures for the clusters in our sample which exhibit significant 
cooling. These temperatures were then employed in the spectral determination of integrated 
luminosities for these objects. Overall, our sample displays temperatures in the range of 
3.4 — 12.1 keV, abundances which range from 0.1 — 0.8 times the solar value, and unabsorbed 
bolometric luminosities within R2500 which span 4.9 — 56.1 x 10^^ erg s~^. 

Cluster gas and total masses within R2500 and R200 were estimated using the outcomes 
of spectral and surface brightness fitting, resulting in virialized cluster masses of 2.7 — 21 x 
10^^ /ifo^ Mq (Table ED and respective cluster gas masses of 0.13 - 2.9 x 10^^ /i^ M©. A 
weighted average of ecas mass fractions e;ives fgas = 0.136 ± 0.004 h-fo , res ulting in a 



0.28=fc 0.01 (68% confidence, formal error), this value is in good agreement with lSpergel et al. 
fl2007h . 



Dynamical masse s within R2nn were ca lculated via the Jeans Equation using velocity 
dispersions taken from lBorgani et al. I (119991 ) (Table [8]). Comparisons between X-ray masses 
and dynamical masses result in a weighted average of Mdyn/Mx-ray = 0.97 ±0.05, indicating 
good agreement between between these two methods, however a fair amount of scatter is 
evident (~ 30%; Figure ED- Dynamical masses are noticeably larger in the case of clusters 
which exhibit significant substructure (MS1006+12 and MS1008-12), a factor which may b e 
responsible for an overestimation of velocity dispersions (IBird. Mushotzky. &: Metzlerlll995l ). 



Weak lensing masses within 500 h^QQ kpc were obtained for seven of the clusters in our 
sample. X-ray masses were calculated within this region and compared to the lensing results 
(Table [9]). Although the distribution appears somehwat asymmetric (Figure [6]), a weighted 
average gives Miens/Mx-ray = 0.99 ± 0.07, with a scatter of ~ 30%. 



X-ray scaling laws fo r this s ample were investigated in a manner similar to lAUen. Schmidt, fc Fabian 



( I2OOII ) and lEttori et al. I (120041 ). taking cosmological factors into account. The best fitting 



Lx — Tx relationship for our sample results in an intercept of 0.13^° and a slope of 



2.4 ± 0.2 (Figure [7]). While this slope is lower than some estimates (lArnaud fc Evrardlll999 
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Ettori et al. Il2004) . it is consistent with both lAUen. Schmidt, fc Fabian! (1200 ll ) and lEttori. De Grandi. fc M 
( 120021 ) ■ Moderate negative evol ution is ind icated by this index being greater than that 
predicted by self-similar models (j VoitI 120041 ) . The Mass-Tx relationship for these clusters 
also exhibits a somewhat steeper slope than expected, at 1.7 ± 0.1, with an intercept of 
l.Stoi X 10^=^ (Figure ED. 



The best fitting scaling laws for our sample (listed above) were combined with Equa- 
tion [ini with the ultimate goal of calibrating relationships between red-sequence optical 
richness (Bgc^red) and global cluster parameters {Tx, Lx and M25oo)- We find that Tx scales 
relatively well to Bgc,red with a 40% scatter (Figure [9]). M2500 shows a scatter of ~ 70% (Fig- 
ure [11]), which is consistent with the Tx-Bgc,red scatter given that M cx T^/^. The Lx-Bgc,red 
relationship exhibits a significantly larger scatter at ~ 200% (Figure [TUl) . 

Our results indicate that Bgc^red does exhibit intial promise as a mass indicator. Accurate 
calibration of a relationship between optical richness and cluster mass will require the use 
of additional clusters which posess both well-constrained X-ray temperatures and optical 
richness measurements. 

Overall we find that multiple cluster mass estimators - dynamics, weak lensing and 
X-ray observations, along with optical richness in the cluster red sequence - are converging 
for this sample of well-studied clusters. While individual correlations still have significant 
scatter, there is little evidence for large systematic bias in any of these methods. Cluster 
characteristics which might be considered problematic for one or more techniques (cluster 
substructure, merging, and/or the presence of a cool core) appear to perturb these relation- 
ships relatively little, as long as high quality data are obtained and the analysis is tuned 
to correct for these. In particular the correlation between cluster temperature and optical 
richness, the most easily obtained of the mass estimators, is promising. A remaining concern 
is that X-ray selected clusters may not prove to be typical of all massive clusters at these red- 
shifts. Additional checks involving clusters covering a broader range of redshift and selection 
technique (e.g., SZ, optical, weak lensing) will be necessary to solidify our understanding of 
the most reliable and efficient methods of cluster mass estimation. 
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Fig. 1. — Optical Images. X-ray contours are overlaid on optical images of clusters in 
our sample. HST images of Abell 2390, MS0016+16, MS0302+16, MS0440.5+02, MS0451- 
03, MS1358+62, MS1455+22, MS1512+36, and MS1621+26 were retrieved from the MAST 
archive, with respective filters F814W, F555W, F814W, F702W, F702W, F814W, F606W, 
F555W, and F814W. Gunn r band (to r ~ 24) CFHT images for the remaining five clusters 
('MS045 1+02, MS0839+ 29, MS0906+11, MS1006+12, and MS1008-12) are from the CNOCl 
survey (lYee et al. Ill996l ). X-ray countours were created from adaptively smoothed Chandra 
0.6-7.0 keV flux images, and all have linear values except for three of the cooling core 
clusters (Abell 2390, MS0440+02, and MS1455+22) which are overlayed with logarithmic 
X-ray contours. Note the significant amount of substructure present in MS0451+02. 
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Fig. 2. — Surface Brightness Profiles. Radial surface brightnesses of CNOC clusters 
calculated in 1" annular bins using 0.6-7.0 keV data. Horizontal dotted lines denote the 
best-fitting background of the model. Pairs of vertical dotted lines indicate regions of the 
profile which were excluded from fitting, while a single vertical dotted line denotes the 
radius to which fitting was performed. Solid lines describe the best fitting single f3 model 
for each cluster, while short dashed lines indicate double f3 model fits to the data. Long 
dashed vertical lines represent R25oo- Cluster emission is well described by either a single or 
double f3 model except in the case of MS0451+02, which exhibits a significant amount of 
substructure. Seemingly high values of reduced are a consequence of fitting high resolution 
data with small error bars, such that even small deviations from symmetry will drive up 
significantly. 
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Fig. 3. — Temperature Profiles. Radial temperature profiles determined by fitting, in 
XSPEC, single temperature spectral models to spectra extracted from circular annuli selected 
to include a minimum of 2500 counts in the 0.6-7.0 keV band. Error bars represent 90% 
confidence limits. Indications of cool cores are seen in Abell 2390, MS0440-I-02, MS0839-I-29, 
MS1358+62, and MS1455+22. 
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Fig. 4. — Universal Temperature Profiles. Top - Individual temperature profiles of the 
five cooling core clusters in this sample. Abell 2390: squares and solid line, MS0440+02: 
diamonds and dashed line, MS 0839+29: filled circles and dot-dash line, MS1358+62: trian- 
gles and dot-dot-dot-dash line, MS455-I-22: x-symbols and long dashed line. The best fit to 
th e data is overlayed as a thick soli d line, while the dotted line indicates the best fit found 
by lAUen. Schmidt. &: FabianI (120011 ) . Bottom - Weighted average temperature profile of the 
combined five clusters in 0.1 R/R2500 bins, with best fitting models overlayed as described 
above. 
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Dynamical vs. X — ray Masses 
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Fig. 5. — Dynamical vs. X-ray Masses. X-ray mass estimates are plotted against the 
ratio MDyn/Mx-ray The dot-dash line indicates a mass ratio of 1.0. Squares represent 
coohng core clusters. Error bars denote 68% confidence limits. 
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Weak Lensing vs. X-ray Masses 
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Fig. 6. — Weak Lensing vs. X-ray Masses. X-ray mass estimates are plotted against 
the ratio MLens/Mx-ray The dot-dash hne represents a mass ratio of 1.0. Squares indicate 
coohng core clusters, and error bars denote 68% confidence limits. Though the distribution 
appears assymmetric, it is not statistically significant. 
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0.4 0.6 0.8 1.0 1.2 

Logio Tx [keV] 



Fig. 7. — The Lx — Tx Relationship. X-ray temperatures are plotted against cosmologi- 
cally corrected unabsorbed bolometric luminosities within R25oo- Squares represent clusters 
with cool cores. The solid line indicates the best fitting relationship, with a slope of 2.4±0.2. 
The dash-dot line illustrates a slope of 2.0. Error bars indicate 90% confidence limits. 
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Fig. 8. — The Mass-Tx Relationship. X-ray temperatures are plotted against X-ray 
derived masses within R25oo- Squares represent clusters with cool cores. The solid line 
indicates the best fitting relationship, with a slope of 1.7 ± 0.1, while the dash-dot line 
represents a slope of 1.5. Error bars indicate 90% confidence limits. 
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Fig. 9. — Tx vs. Bgc^red- A log- log plot of Tx vs. Bgc^red for the clusters in our sample. Error 
bars represent 68% confidence intervals. The solid line indicates the best fitting relationship 
between t hese two paramete r s (Se ction [9l), while the dash-dot line illustrates the previous 
results of lYee and EUingsonl (120031 ). Squares denote clusters which exhibit significant cool 



cores. 
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Fig. 10. — Lx vs. Bgc^red- A log-log plot of Lx VS. Bgc^red for the clusters in our sample. Error 
bars represent 68% confidence intervals. The solid line indicates the b est fitting relationship 
(Section[9]), and the dash-dot line denotes the lYee and EUingsonl (l2003l ) result. Squares again 
denote clusters which exhibit significant cool cores. 
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Fig. 11. — X-ray Mass vs. Bgc,red- A log-log plot of X-ray mass estimates vs. Bgc,red 
for the clusters in our sample. Error bars represent 68% confidence intervals. The solid 
line represents the best fitting relati onship between th e se tw o parameters (Section [9]). The 
dash-dot line was produced from the lYee and EUingsonl (120031 ) results, scaling their intercept 
by the average of M2500 /M200 for this sample. Squares denote coohng core clusters. 
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Table 1. Cluster Sample 





Cluster 




z 


Centroid 




L" 


obsid 


Array 


Chandra Exposure* 










[RA] 


[Dec] 




kpc] 






[seconds] 


Abell 2390 


0, 


.2279 


21:53:36.794 


+17:41:41.85 


3, 


.65 


4193 


ACIS-S 


89624 


MS 


0015.9+1609 


0, 


.5466 


00:18:33.64 


+ 16:26:11.2 


6, 


.39 


520 


ACIS-I 


69235 


MS 


0302.7+1658 


0, 


.4246 


03:05:31.72 


+17:10:01.5 


5, 


.57 


525 


ACIS-I 


11764 


MS 


0440.5+0204 


0, 


.1965 


04:43:09.974 


+02:10:18.01 


3. 


.17 


4196 


ACIS-S 


45104 


MS 


0451.5+0250 


0, 


.2010 






3. 


.31 


4215 


ACIS-I 


66275 


MS 


0451.6-0305 


0, 


.5392 


04:54:11.19 


-03:00:52.2 


6, 


.34 


902 


ACIS-S 


43652 


MS 


0839.8+2938 


0, 


.1928 


08:42:55.999 


+29:27:25.45 


3, 


.21 


2224 


ACIS-S 


31383 


MS 


0906.5+1110 


0, 


.1709 


09:09:12.81 


+10:58:31.7 


2, 


.91 


924 


ACIS-I 


31392 


MS 


1006.0+1202 


0, 


.2605 


10:08:47.56 


+11:47:34.0 


4, 


.03 


925 


ACIS-I 


25590 


MS 


1008.1-1224 


0, 


.3062 


10:10:32.44 


-12:39:41.4 


4. 


.52 


926 


ACIS-I 


37376 


MS 


1358.4+6245 


0, 


.3290 


13:59:50.640 


+62:31:04.20 


4. 


.74 


516 


ACIS-S 


53055 


MS 


1455.0+2232 


0, 


.2570 


14:57:15.110 


+22:20:32.26 


3, 


.99 


4192 


ACIS-I 


91886 


MS 


1512.4+3647 


0, 


.3726 


15:14:22.507 


+36:36:20.15 


5, 


.14 


800 


ACIS-S 


14665 


MS 


1621.5+2640 


0, 


.4274 


16:23:35.37 


+26:34:19.4 


5, 


.59 


546 


ACIS-I 


30062 



Corrected exposure time (see text for corrections applied). 
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Table 2. /?-Model Fits 



Cluster 






-'0 




-'B 


XVDOF 


Abell 2390 


631^ 


n cqi +0.003 
U.J01_o 003 


15751^^ 




1 7+0.2 
' -0.2 


368.8/145 


MS 0015.9+1609 


2531}^ 


n 79+0.02 
' ^-0.02 


lOOl^ 


1 


1 7+0.04 
' -0.04 


316.2/230 


MS 0302.7+1658 


27+^ 


U-JU_o 02 


3371^^ 





oQ+0.04 
•^y-0.04 


328.2/324 


MS 0440.5+0204 




n 444+0-006 

U.+++_o.006 


qqq+84 




2.41°:} 


372.2/191 


MS 0451.5+0250 


5201^° 


0-87to:o8 


511} 




1 7+0.3 
' -0.3 


620.4/227 


MS 0451.6-0305 


200+}° 


n 70+0.02 
"J- ' '-'-0.02 


1541^ 


1, 


01 +0.05 
•^J-0.05 


257.3/205 


MS 0839.8+2938 


58+^ 


0.621°:°? 


4071^^ 


2 


qq+0.07 
■ "5 '-•-0.07 


220.7/208 


MS 0906.5+1110 


96l6^ 




185+? 




1 0+0-2 
-'-■^-0.2 


177/174 


MS 1006.0+1202 


1291}^ 


U.OO_o.o2 






0.81°:? 


270.7/213 


MS 1008.1-1224 


1131^ 






1 


Qq+0.06 

■'-'^-0.06 


329.4/218 


MS 1358.4+6245 


76+5 


fi4+0-0i 
'-'•^+-0.01 


455li 




3.81°:} 


349.6/207 


MS 1455.0+2232 


36.7l°i 


612+" °°'^ 

U.UiZ_o.0O4 


3277+^^ 


2 


90+0.06 
.^^-0.06 


789.1/209 


MS 1512.4+3647 


22+3 


51+0 01 

'-'•^-L-O.Ol 


6881^- 


1 


co+O.OS 
.'-'^-0.08 


192.5/191 


MS 1621.5+2640 


18411 


cifi+0-05 


3411 


1, 


00+0. 09 
."^^-O.l 


234/225 



Surface brightness / in units of 10 ^ photons sec ^ cm ^ arcsec 



Table 3. Double /3-Model Fits 



Cluster 


[/iyo^kpc] 


Pi 




i?c2 [^70^kpc] 


P2 


-'02 


-'B 


XVDOF 


Abell 2390 
Abell 0440.5+0204 
MS 0839.8+2938 
MS 1358.4+6245 
MS 1455.0+2232 


20^^ 

2111 
10+2 
36l« 


"-'•'-''-'-0.07 

i.olSi 

o.9;°i 

'-'•^'-'-0.02 
-'-•-'--0.2 


17091^^6^ 
12461}^^ 

7f5i-H146 
"-"-■--lOS 

2635_383 


331:1 

59+8 
561?^ 


0.521°:°? 
0.4891°:°^} 
62+°-°^ 

'-'•'-'^-0.02 

63+°-^ 
64+°-°^ 


504.9lif 
3381?? 
366+^? 
188+i 

13061^?^ 


-1 7+0.2 
-■-• ' -0.2 

2.80lo:o9 

9 q+0.1 
^•'-'-0.1 

^•"-'-0.4 
9 4cr+0.07 
^•^^-0.07 


177.6/106 
221.5/191 
190/167 
156/139 
276.3/183 



^Surface brightness / in units of 10 ^ photons sec ^ cm ^ arcsec 
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Table 4. Single Temperature Spectral Fits for R < R2500 





Cluster 


R-2500 


kT 


Nh 


Z 


XVDOF 






[hjQ kpc] 


[keV] 


[10^0 cm-2] 


[solar] 




Abell 2390 






6.8 




401.8/434 


MS 


0015.9+1609 


550l| 




4.07 




185.5/271 


MS 


0302.7+1658 




4+3 


10.9 




10.3/26 


MS 


0440.5+0204 


640+^f 




9.67 




272.2/378 






583^^^ 


7+2 
' -2 


121^ 




267.7/377 


MS 


0451.5+0250 


496+?° 


"•-■--0.6 


7.8 




443.9/420 


MS 


0451.6-0305 




10.2+i 


5.07 




223/258 


MS 


0839.8+2938 


4911^^ 


4 n+o-3 


4.11 


U.U_o 2 


241.7/283 


MS 


0906.5+1110 


5821^^ 




3.54 


2+°-2 


134.2/319 


MS 


1006.0+1202 


557_^g 




3.76 


3+°-2 

U.O_Q 2 


100/199 


MS 


1008.1-1224 


501;^^ 


p: q+0.9 


6.98 


U.Z_Q 2 


159.5/203 


MS 


1358.4+6245 


575+^^ 




1.93 


U.O_Q 2 


272.2/332 


MS 


1455.0+2232 


496t^ 


4.4;°:} 


3.13 


0.38+n^ 


569.6/385 


MS 


1512.4+3647 


351+i 


3.4+2;? 


1.38 


0.6l°i 


43.2/74 


MS 


1621.5+2640 


493ir 


7tl 


3.57 




56.3/116 
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Table 5. Cooling Core Corrected Single Temperature Spectral Fits for R < R2500 



Cluster 


R-cutoff 

[Ko kpc] 


^2500 

[Ko kpc] 


kT 
[keV] 


Nh 
[10^0 cm-2] 


Z 

[solar] 


XVDOF 


Abell 2390 
MS 0839.8+2938 
MS 1358.4+6245 
MS 1455.0+2232 


66 
109 
154 
211 


7611^^ 

658li 
575li 


-'-^•-'--0.9 

9+2 

5.91^:? 


6.8 
4.11 
1.93 
3.13 


3+° -'^ 
6+°-^ 

U.U_Q 4 

0.4lg;i 

0.41^:^ 


355/434 
195.5/237 
209/253 
238/273 
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Table 6. Mass Estimates for R < R2500 



Cluster 


no 


Mgas 


M2500 




[10-1 cm-3] 




[10^3 hfo^ Mq] 


Abell 2390 


n 224+°-°°^ 


e.ilSi 


591^ 


MS 0015.9+1609 


0.0958+o;ooo8 


4.61°:} 


38I3 


MS 0302.7+1658 


0.461H1 


0.81°:} 


q+3 


MS 0440.5+0204 


0.2801°:°°! 


1 62+°°^ 


26+^ 


MS 0451.5+0250 


0.03411°:^! 


9 0+0.3 


17.61°:^ 


MS 0451.6-0305 


n 1Q9c;+0.0009 
U.lOZO_Q Q009 


c 9+0.2 
^•^-0.1 


521! 


MS 0839.8+2938 


n 24n+°-°°^ 


1 6+°-i 

-■-•"-0.1 


231^ 


MS 0906.5+1110 


n 190+0.001 

U.iZO_o 001 


2 00+° °^ 


24I2 


MS 1006.0+1202 


n nyq+o.ooi 
u.u^y_o 001 


9 0+0.1 
^."-•-0.1 


25li 


MS 1008.1-1224 


100+°°°^ 


2 0+°-i 
^••J-o.i 


191^ 


MS 1358.4+6245 


0.1391°:°°^ 


4+0.5 


471^ 


MS 1455.0+2232 


0.501°:°} 


9 7+0.3 
^. ' -0.3 


271^ 


MS 1512.4+3647 


6S+°-°i 


0.861°:°^ 


71{ 


MS 1621.5+2640 


A Qcro+O.OOl 




241^ 
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Table 7. Mass Estimates for R < R200 



Cluster 


a 
rv200 




M200 


f 

J QCLS 




[hfo' Mpc] 


[10^^ h,-'/^ Mo] 


[10^^ hfo' Me] 




Abell 2390 


9 7+0.1 
^- ' -0.1 


q+0.3 

"-"•^-O-S 




0.191°:°^ 


MS 0015.9+1609 


9 1+0.1 
^•J--o.i 


9 «9+0.l0 
^•'J^-o.io 




1 s+°-°^ 

•J- -'-'-'-0.01 


MS 0302.7+1658 


1 9+0.5 
-^•^-0.3 


56+°-°'^ 




17+006 

-0.06 


MS 0440.5+0204 


9 1+0.3 


1 14+0-07 




1 3+°°2 

'-'■J^'J-0.02 


MS 0451.5+0250 


9 p;+0.2 
^•'-'-0.1 


9 0+0.3 


i5.9i°:« 


1 7+°°2 
'-'■-'-'-0.02 


MS 0451.6-0305 


9 q+0.2 
^•"-"-O-l 


9 qq+0.10 
^•'J'-'-O.lO 




'-'•-'--L-O.Ol 


MS 0839.8+2938 


2 0+°-3 

^•U-0.2 


n 70+0.06 

^- ' "J-O-OO 




oq+°°^ 


MS 0906.5+1110 


9 1+0.2 
^•-■--0.1 


J--^^-0.05 


8.61°:? 


i2+°-°i 


MS 1006.0+1202 


2.0l°J 


1 6+°-^ 


91} 


^•-i^ ' -0.02 


MS 1008.1-1224 




-'-•'J^-0.07 


7 9+0.6 

' -^-0.6 


0.191°:°^ 


MS 1358.4+6245 


9 4+0.3 


1 7+0.3 
-■-•'-0.3 


171^ 


0.091°:°^ 


MS 1455.0+2232 


2.04+°? 


1 9+0.1 
-■-•^-0.1 


9.51^:^ 


1 2+°°^ 
'-'■-'-'^-0.02 


MS 1512.4+3647 


1.31°:1 


0.591°:°^ 


9 7+0.3 


A 99+0.03 
'-'■■^'^-0.03 


MS 1621.5+2640 


1-81^:^ 




101^ 


'-'■^^-0.03 



^90% confidence intervals 
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Table 8. Dynamical Mass Comparisons for R < R200 



Cluster 


Dynamical Mass 


X-ray Mass 


Ratio 




[10^^ hfo^ Mq] 


[10^^ hfo^ Mo] 


[Dynamical / X-ray] 


Abell 2390 


261^ 


2lli 


-'-■^-0.2 


MS 0015.9+1609 




17+1 


0.9l°i 


MS 0302.7+1658 




3li 


1 q+0.5 


MS 0440.5+0204 




91^ 


0.7lgJ 


MS 0451.5+0250 


151^ 


15.9l°i 


0.91°;} 


MS 0451.6-0305 


241^ 


211^ 


l-ll^J 


MS 0839.8+2938 


12l3 


81} 


1-51^:1 


MS 0906.5+1110 


81? 


8.61°;? 


0.91°;? 


MS 1006.0+1202 


13+3 


91} 


1 4+0-5 

-■-■^-0.4 


MS 1008.1-1224 


12+3 


7 r,+0.6 
' -^-0.6 


1-71^:^ 


MS 1358.4+6245 


141? 


171^ 


0.81°;? 


MS 1455.0+2232 


12li 


9.5l°| 


1 q+0.3 

^■"^-0.2 


MS 1512.4+3647 


2-61J.8 


9 7+0.3 
^- ' -0.4 


i.olSi 


MS 1621.5+2640 


811 


101^ 


0.81°;^ 



Table 9. Weak Lensing Mass Comparisons for R < 500 h ^ kpc 



Cluster Weak Lensing Mass X-ray Mass Ratio 

[10^^ h-i Mq] [10^^ h-i Mq] [Lensing/X-ray] 



Abell 2390 




c; q+0.5 
0.0_Q 4 


0.9+°;} 


MS 0015.9+1609 




c 9+0.4 
^■^-0.4 




MS 0906.5+1110 


o c+0.7 
0.0_Q 7 


9 O+0.2 


1 9+0-3 
^■^-0.3 


MS 1358.4+6245 


4 2+0-8 

^■^-0.8 


4 q+0.6 

^■^-0.6 




MS 1455.0+2232 


q 9+0.7 
"^■^-0.7 


q 9+0.3 
^■■^-0.3 


1 n+°-2 

-'-•U-0.2 


MS 1512.4+3647 


2.ol°:« 




1 4+0-6 
-■-•^-0.6 


MS 1621.5+2640 


51} 


3.61^:? 


1 4+0.4 
-■-■^-0.4 
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Table 10. Cluster Richness and Luminosity 





Cluster 


[10^^ erg s"^] 


[Mpc^-^^] 


Abell 2390 


56.ll°:f 


1473 + 218 


MS 


0015.9+1609 


41.0+°-^ 


1789 ± 260 


MS 


0302.7+1658 




762 ± 201 


MS 


0440.5+0204 




559 ± 154 


MS 


0451.5+0250 




1198 + 208 


MS 


0451.6-0305 


49.i;r 


1232 + 180 


MS 


0839.8+2938 


r 07+0.1 b 
^•^ ' -0.08 


1106 + 187 


MS 


0906.5+1110 


o-'-'-O.S 


1134 + 209 


MS 


1006.0+1202 




^•^-0.6 


1330 ± 216 


MS 


1008.1-1224 


9.0l°j 


1444 ± 212 


MS 


1358.4+6245 


17.7l°f 


1182 + 205 


MS 


1455.0+2232 


90 o+O lb 
■^'J--^_0.2 


518 + 154 


MS 


1512.4+3647 


4.9l°i 


582 ± 180 


MS 


1621.5+2640 




996 ± 199 



'^Unabsorbed bolometric X-ray luminosity for 
R < R2500 

''Cooling core corrected 
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Table 11. Fitting Parameters (A = 2500) 



Relationship 




C2 


Tx — Bg(..red 


-2.4 ±0.8 


1.1 ±0.3 


E{z) Lx — Bgcred 


-6 ±2 


2.3 ±0.7 


E{z) M2500 — Bgc,red 


9 ±2 


1.8 ±0.5 



